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ABSTRACT 


Mass  resolved  excitation  spectra  are  presented  for  the  benzyl  radical  and  its  clusters  with 
Ar,  N2,  CH4,  C2H6,  and  C3Hg.  The  cluster  spectra  exhibit  small  red  shifts  (<  50  cm'1)  relative  to 
the  unclustered  benzyl  radical  for  the  Dj(12A2)  4-  Dq(12B2)  and  D2(22B2)  4—  DQ(12B2)  electronic 

transition  regions.  A  unique  set  of  low  energy  van  der  Waals  modes  is  observed  for  these  clusters 
for  each  excited  electronic  state  investigated.  The  cluster  spectra  also  reveal  significant  vibronic 
coupling  between  the  two  excited  electronic  states  of  the  benzyl  radical,  as  evidenced  by  a  single 
vibrational  predissociation  threshold  for  each  cluster.  Ab  initio  calculations  are  performed  on  the 
benzyl  radical  to  examine  excited  electronic  state  structure,  predict  transition  energies,  estimate 
ionization  energy,  and  determine  partial  atomic  charges  in  the  electronic  states  of  interest.  The 
resulting  partial  charges  are  used  in  empirical  atom-atom  potential  energy  calculations  to  aid  in  the 
understanding  of  cluster  spectroscopic  shifts,  binding  energies,  and  van  der  Waals  modes. 
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INTRODUCTION 


Most  of  the  spectroscopic  information  on  reactive  intermediates  accrues  from  either  matrix 
isolation  or  room  temperature  gas  phase  experimental  studies.  Within  the  last  ten  years, 
however,  studies  of  reaction  intermediates,  in  particular  radicals,  have  been  reported  employing 
supersonic  jet  expansion  techniques.  These  recent  radical  generation  techniques  involve  either 
electric  discharge,  photolysis,  or  thermal  decomposition  of  suitable  precursors,  occurring  as  part  of 
the  supersonic  expansion.  When  such  techniques  are  combined  with  fluorescence,  mass  resolved 
excitation,  or  dispersed  fluorescence  detection,  very  detailed  information  can  be  obtained 
concerning  the  ground  and  excited  state  properties  of  radicals. 

Over  the  past  35  years  the  benzyl  radical  has  been  the  subject  of  many  experimental  and 
theoretical  studies.24'32  The  benzyl  radical  is  a  simple  aromatic  stabilized  7t  radical  with  relatively 
high  symmetry  (C2v).  Early  theoretical  calculations  correctly  predict  the  ground  state  to  have  B2 

symmetry,  but  disagree  as  to  the  ordering  of  the  1  A2  and  2  B2  first  two  excited  states.  This 

theoretical  uncertainty  as  to  the  ordering  of  the  excited  states  is  understandable  considering  their 
close  proximity.  (It  is  now  known  that  the  22B2  state  lies  roughly  400  cm'1,  or  0.05  e  V,  above 
the  1  ^  state.)  Experimental  work  conducted  by  Leach  and  co-workers  reveals  that  a  low  lying 
first  excited  (Dj)  state  of  A2  symmetry  lies  at  22,000  cm"  .  Their  determination  of  the  D1  excited 

state  symmetry  relies  on  a  rotational  band  contour  analysis  of  the  emission  spectra  of  several 
deuterated  benzyl  radical  species.  This  assignment  was  later  confirmed  by  Friedricht  et.al.  An 

inspection  of  the  benzyl  radical  absorption  and  emission  data  reveals  the  existence  of  three  intense 
features  between  300-650  cm  above  the  assigned  1  A2  electronic  origin  with  no  readily 

identifiable  comparable  features  accessible  in  the  ground  state.  Cossart-Magos  and  Leach19 

performed  a  linear  pseudo-Jahn-Teller  (LPJT)  vibronic  coupling  calculation  which  treats  the 
interaction  between  the  1 8b  and  6b  vibrations  of  the  1A2  electronic  state  with  the  2  B2 

electronic  state.  Their  calculations,  which  account  for  the  experimentally  observed  isotope  effects, 
predict  the  unperturbed  22B2  electronic  origin  to  lie  roughly  400  cm'1  above  the  12A2  origin. 

More  refined  vibronic  coupling  calculations  have  been  performed  ’  which  include  the  8b  ,  in 


addition  to  the  18b1  and  6b1  vibrations.  These  reanalyses  of  the  problem  predict  the  unperturbed 
2  B2  electronic  state  to  be  -850  cm  above  the  1  A2  electronic  state.  Recently,  well  resolved 

supersonic  jet  generated  spectra  of  the  benzyl  radical  have  been  obtained  using  time-of-flight  mass 
resolved,9  and  dispersed  emission11  (DE)  experimental  techniques  which  confirm  the  location  of 
these  vibronic  features. 

Although  radicals  are  known  to  play  an  important  role  as  reaction  intermediates  in 
condensed  phases,  little  is  understood  about  their  solvation  physics  and  chemistry.  Here  we 
report  a  study  of  the  effect  of  solvation  by  non-polar  molecules  on  the  electronic  transitions  of  the 
benzyl  radical.  Many  of  the  phenomena  observed  for  the  radical  solute-solvent  clusters  examined 
here  (e.g.,  small  spectroscopic  red  shifts,  appearance  of  van  der  Waals  vibrational  modes  in  the 
cluster  spectra,  vibrational  predissociation,  vibronic  coupling)  are  known  to  be  important  for 
clusters  of  stable  species. 

The  experimental  results  presented  in  this  report  include  cluster  transition  energy  shifts  with 
respect  to  the  bare  radical,  an  estimation  of  cluster  binding  energy  for  different  radical  electronic 
states,  and  an  identification  of  cluster  vibrational/torsional  van  der  Waals  structure.  Calculations  of 
benzyl  radical  electronic  states,  cluster  structure,  and  binding  energies  are  employed  to  aid  in  the 
elucidation  of  the  experiment  results.  In  particular,  ab  initio  calculations  are  performed  to  examine 
electronic  state  structure  and  to  determine  atomic  partial  charges  in  the  ground  and  first  two  excited 
electronic  states.  These  partial  charges  are  used  in  conjunction  with  dispersion  and  repulsion  atom- 
atom  interaction  parameters  to  estimate  cluster  structure,  binding  energy,  and '  ibrational/torsional 
van  der  Waals  motion. 

II.  PROCEDURES 

1  Experiments 

Experiments  based  on  the  generation  and  detection  of  radicals  are  similar  to  those  of  stable 
species  performed  previously  in  this  laboratory,34  although  subtle  differences  are  encountered  in 


studying  the  two  types  of  samples.  What  follows  is  a  discussion  of  these  differences  as  well  as  a 
detailed  description  of  the  formation  and  detection  of  the  benzyl  radical  and  its  clusters. 

Figures  1  and  2  pictorially  represent  the  optical  and  electronics  configuration  used  in  our 
experimental  study  of  radicals.  Figure  1  schematically  illustrates  the  generation  and  detection  of 
radicals.  An  ArF  excimer  laser  with  an  output  power  of  roughly  80  mJ/pulse  was  focused  onto  the 
exit  region  of  an  R.M  Jordan  pulsed  valve  operating  at  10  Hz.  Helium  was  used  as  the  expansion 
gas  at  a  backing  pressure  of  80  pisg  and  clusters  were  formed  by  using  an  expansion  mixture 
containing  10%  of  the  clustering  agent  of  interest.  The  radicals  and  their  clusters  generated  in  the 
high  density  photolysis  region  are  collisionally  cooled  and  travel  downstream  to  a  pair  of  ion 
extraction  plates  (flight  time  roughly  50|is)  where  they  were  optically  excited  to  a  particular 
vibronic  state  of  Dj  or  D2  and  ionized  by  the  absorption  of  a  second  photon.  The  ions  created 

were  extracted  by  a  350  V  potential  difference  and  accelerated  into  a  time-of-flight  mass 
spectrometer  by  -4  keV.  Mass  resolved  detection  of  the  ionized  radical  species  was  accomplished 
via  a  Galieo  Electro-Optics  multi-channel-plate  (MCP)  detector  at  the  end  of  a  1.5  m  flight  tube, 
generating  mass  resolved  excitation  spectra  (MRES).  For  the  systems  presently  under  study,  two 
Nd:YAG  pumped  dye  lasers  were  used  as  the  excitation  and  ionization  source.  Either  C460  or 
C440  was  used  as  the  dye  lasing  medium  for  the  Dj,D2  <—  DQ  excitation  source  and  F548 

(frequency  doubled)  was  used  as  the  dye  lasing  medium  for  the  I «-  Dj,D2  ionization  source. 

The  configuration  of  electronic  instruments  needed  to  control  the  timing  of  the  above  events 
is  schematically  illustrated  in  Figure  2.  A  trigger  pulse  originating  from  the  oscillator  of  laser  1 
was  sent  to  the  oscillator  of  laser  2  such  that  the  Q-switches  of  lasers  1  and  2  were  activated 
simultaneously.  This  "common"  oscillator  trigger  was  also  passed  through  a  delay  stage  and  used 
to  trigger  the  pulsed  valve  power  supply.  The  output  trigger  from  the  pulsed  valve  power  supply 
was  then  sent  to  another  delay  stage  and  used  to  activate  the  ArF  excimer  laser.  Detection  of  the 
radicals  was  achieved  by  sending  the  variable  synch  trigger  of  laser  1  to  a  Stanford-Research- 
System  (SRS)  delay  box  which  controls  both  the  boxcar  averager/oscilloscope  timing  and  the  data 
acquisition  peripherals.  Since  the  firing  of  the  excimer  laser  created  a  0.9  V  RF  noise  pulse 


transmitted  through  our  detection  electronics,  the  input  threshold  on  the  SRS  instrument  was  set  at 
greater  than  1.0  V.  To  alleviate  additional  grounding  and  feedback  problems,  an  opto-isolator  was 
placed  between  the  Q-switch  of  laser  2  and  the  electronics  of  the  wavelength  extension  system 
(WEX)  in  order  to  amplify  the  triggering  pulse. 

The  next  step  taken  in  achieving  benzyl  radical  signal  was  to  align  (spatially  and  temporally)  the 
excitation  laser,  ionization  laser,  and  photolysis  laser.  This  is  accomplished  through  a  four  step 
procedure.  First,  the  excimer  laser  output  (193  nm)  was  spatially  aligned  by  removing  the  back 
plate  of  the  vacuum  chamber  and  focusing  the  light  beam  to  a  spot  size  of  5  mm  onto  the  exit 
region  of  the  nozzle.  Second,  the  output  from  the  excitation  and  ionization  lasers  was  reduced  in 
intensity  and  focused  onto  a  fast  photodiode  (EG&G  FND100  reversed  biased  by  15  V,  response 
time  <5ns)  with  the  signal  output  sent  to  an  oscilloscope.  This  allowed  a  temporal  overlap  of  the 
two  laser  pulses  to  be  achieved.  Third,  a  one  color  mass  detected  aniline  signal  (at  roughly  36,400 
cm'1)  was  used  to  optimize  the  nozzle  delay  and  spatial  alignment  of  the  ionization  laser.  The 
aniline  signal  was  also  used  as  a  check  of  mass  calibration.  Finally,  the  excimer  laser  was 
temporally  aligned  while  monitoring  the  one-color  aniline  signal:  at  a  unique  time  delay  the  firing 
of  the  excimer  resulted  in  a  decrease  of  the  observed  aniline  signal  by  more  than  a  factor  of  10. 
This  signal  decrease  is  attributed  to  photolysis  of  aniline  (possibly  multiphoton)  into  unknown 
products.  Since  the  only  remaining  degree  of  freedom  in  the  experimental  set-up  was  the  spatial 
alignment  of  the  excitation  laser,  the  benzyl  radical  signal  could  be  readily  found.  During  optimum 
operation,  2  V  of  benzyl  radical  signal  was  obtained,  with  cluster  signals  varying  from  30  mV  to 
400  mV.  Benzyl  chloride  (stated  purity  99%),  purchased  from  Aldrich  Chemical  Company,  was 
used  as  the  benzyl  radical  precursor.  Negative  peaks  in  the  observed  spectra  depicted  in  the  figures 
result  from  detector  overload  due  to  intense  signals  in  lower  mass  channels. 

2.  Ab  initio  Calculations 

In  order  to  analyze  and  understand  the  benzyl  radical  cluster  spectra  discussed  here  we 
perform  ab  initio  calculations  to  explore  the  effect  of  electronic  excitation  on  the  benzyl  radical.  We 


are  particularly  interested  in  determining  the  electronic  charge  distribution  in  the  ground  and  first 
two  excited  electronic  states.  Calculations  are  carried  out  using  the  HONDO  molecular  orbital 
package  running  on  an  IBM  risc/6000  computer  system  equipped  with  1.5  GB  of  storage  space 
and  64  MB  of  RAM.  Complete-active-space-self-consistent-field  (CASSCF)  calculations  using  the 
Dunning-Huzinaga  [9s,5p]/[3s/2p]  double  zeta-valence  basis  set  (DZV)  are  performed  on  the 

ground  and  first  two  excited  electronic  states  of  the  benzyl  radical.  The  active  space  consists  of 
benzyl  radical  jc  orbitals:  this  restriction  results  in  404  electronic  configurations  for  states  of  B2 

symmetry  and  380  configurations  for  states  of  A2  symmetry.  A  geometry  optimization  was 
performed  for  each  of  the  elctionic  states  by  reducing  energy  gradients  to  less  than  10-4  a.u. 

A 

The  calculations  presented  here  are  similar  to  calculations  performed  by  Rice,  et  al.  who 
chose  the  same  calculational  approach  but  different  basis  sets:  a  full  double  zeta  (DZ)  basis,  and  a 
full  double  zeta  basis  plus  polarization  and  diffuse  functions  on  the  carbon  atoms  (DZd+).  Their 
DZ  results  for  the  three  lowest  electronic  states  of  the  benzyl  radical  are  lower  in  energy  than  our 

DZV  results  by  less  than  0.0005  a.u.  The  DZd+  basis  set  results  give  rise  to  only  a  slight 
improvement  in  12A2*-12B2  transistion  energy.  Since  the  addition  of  polarization  functions  is 

shown  not  to  increase  the  accuracy  of  calculated  partial  atomic  charges  for  small  molecules 
comprised  only  of  carbon  and  hydrogen,  such  as  CH4,  C2H 2,  and  C2H4,  ’  we  have  opted  to 

use  the  DZV  results  to  determine  benzyl  radical  partial  charges.  The  general  procedure  for 

'JO 

obtaining  partial  atomic  charges  has  recently  been  extended  to  larger  systems  and  excited  states. 
What  follows  is  a  brief  description  of  our  method  of  obtaining  partial  charges. 

The  general  method  used  to  estimate  partial  charges  for  a  molecule  involves  creating  a  3- 
dimensional  map  of  the  electrostatic  potential  for  the  electronic  state  of  interest  and  numerically 
fitting  charges  centered  at  the  atomic  positions  to  this  potential.  The  particular  method  we  use  is 
similar  to  that  of  Williams  et  al.36  First,  the  electrostatic  potential  is  evaluated  at  specific  grid 
points.  The  grid  chosen  is  a  simple  3-dimensional  cubic  grid  with  a  spacing  of  2.0  bohr  between 
points.  The  total  grid  size  is  made  sufficiently  large  to  encompass  the  entire  molecule  and  extend 
well  beyond  the  van  der  Waals  radius  of  the  molecule.  Second,  the  electrostatic  potential  points 


used  to  determine  the  partial  charges  are  chosen  by  requiring  that  they  lie  outside  the  van  der  Waals 

radii  of  ail  the  atoms  and  within  a  certain  "shell  thickness".  The  van  der  Waals  radii  used  are39,40 
Rvdw(Carbon)=3.884  bohr  and  Rvdw(Hydrogen)=3.184  bohr.  A  shell  thickness  of  4.0  bohr  is 

employed  for  these  calculations.36  This  value  is  chosen  to  maximize  the  number  of  points  for  the 
determination  of  charges,  while  simultaneously  minimizing  the  errors  introduced  by  assuming  each 
point  has  equal  weight.  Third,  a  standard  non-linear  fitting  routine  is  employed  to  obtain  the  final 
values  of  the  charges41  by  setting  the  magnitude  of  the  error  in  each  point  equal  to  its  value. 

Concern  has  been  expressed  in  the  literature  about  the  type  of  grid  chosen  (i.e.,  cubic, 
elliptical,  etc.).  For  example,  when  benzene's  partial  charges  are  determined  using  a  cubic  grid, 
a  large  discrepancy  between  the  charges  on  carbon  atoms  and  hydrogen  atoms  is  found:  a  4  to  2 
"split"  in  the  carbon  (and  hydrogen)  charges  is  observed.38  While  we  observe  similar  results 
using  various  grid  spacings  and  shell  thicknesses,  our  parameters  above  for  benzene  yield  a 
minimal  4  to  2  split  for  the  carbon  partial  atomic  charges:  carbons  2,  3,  5,  and  6  have  partial  atomic 
charges  of  -0.160  electrons  and  carbons  1  and  4  have  partial  atomic  charges  of  -0.162  electrons. 
The  hydrogen  charges  are  such  that  the  individual  C-H  fragments  have  zero  overall  charge.  For 
methane,  we  calculate  hydrogen  partial  charges  of  +0.138  electrons,  consistent  with  earlier 
studies.  While  a  more  efficient  method  of  generating  electrostatic  grid  points  is  advantageous, 
the  straightforward  approach  employed  here  has  been  shown  to  accurately  predict  electric  multipole 

'in 

properties. 

3 .  Empirical  Atom-Atom  Potential  Energy  Calculations 

Atomic  partial  charges  for  the  benzyl  radical  in  the  ground  and  excited  states  can  be  used  in 

conjunction  with  van  der  Waals  parameters  to  estimate  benzyl  radical  binding  energy  with  the 
solvents  Ar,  N2,  CH4,  C2H6,  and  C3Hg.  Non-bonded  intramolecular  interactions  can  be 

approximated  by  an  atom-atom  pair  potential  function  consisting  of  a  van  der  Waals  exponential-6 
potential  term  and  an  electrostatic  partial  charge  interaction  such  as40 

VBE  =  Z  ij>i  v,j  =  Aifxp(-B||Rij)  •  q/R,/  +  116200  QjQj/Rjj  =  Vvdw  +  V£ 


(1) 


in  which  Ry  in  A  is  the  separation  distance  between  atoms  i  and  j  on  different  molecules,  Ay,  By, 
and  Cy  are  van  der  Waals  parameters,  and  Qit  Qj  in  units  of  electron  charge  are  the  atomic  partial 

charges.  The  numerical  constant  in  Eq.  (1)  is  a  unit  conversion  factor  to  wavenumbers.  Table  I 

lists  values  of  the  van  der  Waals  parameters  used  in  this  study.  For  interactions  between  different 
atomic  species,  and  Cy  are  obtained  from  the  geometric  mean  of  the  homonuclear  interaction 

terms,  whereas  By  is  obtained  from  the  arithmetic  mean.  The  geometry  of  the  benzyl  radical  is 
chosen  as  the  1  B2  CASSCF  optimized  structure,  and  experimentally  determined  geometries  are 
used  for  the  solvents.42  All  solvent  hydrogen  partial  charges  are  taken  as  0.14  e"  with  negative 
charges  equally  distributed  among  the  carbon  atoms  so  as  to  generate  a  zero  net  molecular 

IQ 

charge.  Charges  on  the  nitrogen  molecule  are  placed  on  the  nitrogen  atom  sites  (+0.377  e')  and 
0.25  A  away  from  the  end  of  each  nitrogen  atom  (i.e.,  lone-pair  charges)  on  the  molecular  axis 
(-0.377  e').43  This  distribution  of  charge  reproduces  the  known  nitrogen  quadrupole  moment.44 
Cluster  energies  are  minimized  assuming  rigid  solvent  and  solute  species  using  the  SIMPLEX 
technique  41 

III.  RESULTS  AND  DISCUSSION 

1 .  Benzyl  Radical 

The  2-color  MRES  of  the  benzyl  radical  is  presented  in  Figure  3.9  The  spectrum  contains 
an  origin  for  the  D}  «—  DQ  (12A2  «—  12B2)  transition  at  21,996  cm'1,  intense  features  300-650 

1  2  2 

cm  to  the  blue  of  this  origin  arising  in  pan  from  1  A2  —  2  B2  vibronic  coupling,  and  a  rather 

dense  collection  of  transitions  to  higher  energy.  Summarizing  the  analysis  of  Cossart-Magos  and 
Leach,19  the  two  features  labeled  Aj  and  A3  correspond  to  vibronic  features  of  roughly  an  equal 
hybrid  of  the  22B2  electronic  state  origin  and  the  6b1  vibration  of  the  12A2  electronic  state.  The 
feature  labeled  A2  is  believed  to  be  the  18b  vibration  of  the  1A,  electronic  state  which  remains 

essentially  unchanged  in  energy,  but  borrows  intensity  from  the  2  B2  electronic  state.  The  labeled 
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feature  of  least  intensity  in  this  region  is  the  totally  symmetric  6a  vibration  built  on  the  1  A2 


electronic  state  origin  (Oq).  We  will  comment  further  on  these  assignments  after  analyzing  spectra 
of  the  benzyl  radical  clusters. 

Table  II  contains  results  of  ab  initio  calculations  for  the  benzyl  radical.  The  ground  and 
first  two  excited  states  are  calculated.  A  CASSCF  calculation  using  the  Dunning-Huzinaga  [9s, 
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5p]/[3s,2p]  DZV  basis  set  predicts  the  incorrect  excited  electronic  state  1  A2 —  2  B2  ordering  by 
roughly  1000  cm'1.  The  transition  energies  to  these  states  are  accurately  estimated  to  within 
roughly  1000  cm'1  of  those  observed.  Table  II  shows  that  geometry  relaxation  lowers  the  excited 
state  energies  by  0.07  a.u.,  but  has  little  effect  on  the  ground  state  energy.  The  active  space  for 
these  calculations  consisted  of  the  seven  valence  7t -orbitals  of  the  aromatic  system.  The  major 
electronic  configurations  contributing  to  the  three  states  of  interest  are  found  to  be:  1  B2  state  - 
(b2)2(b2)2(a2)2(b2)1(a2)(b2)(b2)  50%;  Instate  -  (b/  (b/  (a/  (a^1  (b2)  (b2)  (b2)  31%, 
(b2)2  (b2)2  (a^1  (a2)  (b2)2  (b2)  (b2)  22%;  and  22B2  state  -  (b/  (b^1  (a/  (b2)2  (a2)  (b2)  (b2) 
23%,  and  (b2)2  (b2)2  (a2)2  (b2)  (a^  (b2) 1  (b2)  29%.  Other  electronic  configurations  (for  each 
electronic  state)  contribute  less  than  13%  to  each  the  of  the  states.  Also  listed  in  the  bottom  portion 
of  Table  II  is  the  calculated  ionization  energy.  The  calculated  value  of  59,210  cm'1  is  in  reasonable 
agreement  with  our  measured  value  of  58,358  cm'1.9  In  spite  of  the  calculation's  inability  to 
predict  the  correct  state  ordering,  accurate  estimation  of  transition  and  ionization  energies  is 
encouraging. 

As  described  in  the  last  section,  ab  initio  calculations  can  be  used  to  generate  atomic  partial 
charges.  Accurate  atomic  partial  charges  are  necessary  when  evaluating  the  intermolecular  van  der 
Waals  potential  (e.g.,  Lennard-Jones-Coulomb,  and  exponential-6-Coulomb  potentials).  Table  III 
lists  the  results  of  the  partial  charge  calculations  for  the  benzyl  radical  in  its  1  B2,  1  A2,  and  2  B2 

electronic  states.  Listed  in  the  final  column  of  this  table  is  the  calculated  residual  charge,  which  is, 
in  some  sense,  an  indication  of  the  quality  of  the  fit  Note  that  the  a-carbon  partial  charge  is  nearly 
constant  in  all  three  electronic  states  whereas  the  partial  charges  on  carbons  2  and  3  change  by 
more  than  a  factor  of  three  between  the  ground  and  excited  states. 


2 .  Benzyl  Radical(  Ar)  j 

Figure  4  shows  the  MRES  of  the  benzyl  radical(Ar)j  cluster.  Transition  energies  and 

spectral  shifts  for  this  cluster  relative  to  the  bare  benzyl  radical  are  listed  in  Table  IV.  Spectral 
shifts  for  this  cluster  suggest  that  the  feature  labeled  6aQ  is  built  on  the  1  Aj  origin  and  that 

features  A}  and  A2  are  of  a  different  nature.  In  addition,  the  van  der  Waals  vibrational  structure 
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following  A]  and  A2,  which  is  much  less  intense  for  the  1  A2  Oq  and  6aQ  features,  also  supports 
this  assignment.  As  discussed  above,  Aj  and  A2  arise  from  vibronic  interaction  between  the  1  A2 

2  IQ 

and  2  B0  electronic  states,  consistent  with  these  observations. 
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Note  also  that  the  benzyl  radical(Ar)j  spectrum  terminates  at  ca.  450  cm  above  the  1  A2 
origin  as  can  be  seen  by  the  presence  of  the  6aQ  feature  but  not  of  A3.  If  intracluster  energy 

redistribution  (IVR)  and  vibrational  predissociation  (VP)  of  the  cluster  a’-e  fast  on  the  time  scale  of 

our  experiment  (roughly  10  ns),  this  intensity  pattern  suggests  that  the  binding  energy  for  benzyl 
radical(Ar)j  is  greater  than  440  cm'1  and  less  than  650  cm"1.  Furthermore,  the  12A2  and  22B2 

states  must  be  significantly  coupled.  If  this  latter  electronic  state  mixing  did  not  exist,  a  vibronic 
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cluster  spectrum  built  on  the  2  B2  Oq  state  would  extend  for  another  ca.  450  cm  to  the  blue. 

This  spectrum  is  not  observed. 

Exponential-6  potential  energy  calculations  are  performed  for  the  benzyl  radical(Ar)j  cluster 
based  on  the  potential  energy  function  of  eq.(l)  and  the  parameters  given  in  Table  I  for  the  DQ 
1  B2  state.  The  results  are  presented  in  Table  V.  In  order  to  learn  more  about  the  nteraction 

between  the  benzyl  radical  and  the  various  solvents,  we  attempt  to  model  the  excited  state 
radical/solvent  interaction  to  predict  spectroscopic  shifts  and  eventually  cT  itc’-  van  der  Waals 
modes.  The  benzyl  radical  excited  states  are  fit  to  this  potential  form  by  artificially  adjusting  the 
carbon  van  der  Waals  parameters  of  the  benzyl  radical  so  as  to  reproduce  the  benzyl  radical(Ar)j 

spectroscopic  shifts  (see  Table  IV;  this  assumes  binding  energy  is  directly  related  to  transition 
energy  of  the  chromophore).  Specifically,  the  carbon-carbon  interaction  well-depth  for  the  1  A2 

electronic  state,  which  is  proportional  to  the  constants  A  and  C  of  eq.  (1),  is  increased  by  15.3%, 


and  the  well  depth  for  the  2  B2  electronic  state  is  increased  by  8.5%.  This  artificial  adjustment  of 

the  carbon  van  der  Waals  parameters  has  little,  if  any,  physical  justification;  nonetheless,  it  proves 

to  be  instructive  to  examine  if  these  simple  changes  in  parameters  for  the  benzyl  radical  can  roughly 

account  for  the  observed  spectroscopic  shifts,  binding  energies,  and  van  der  Waals 

vibrational/torsional  modes  of  the  remaining  cluster  systems. 

Excited  state  binding  energy  for  the  argon  cluster,  as  listed  in  Table  V,  is  ca.  420  cm'1. 

This  value  appears  to  be  somewhat  low,  since  it  is  outside  the  bracketed  region  of  440-650  cm'1. 

If,  however,  dynamics  of  the  cluster  are  not  sufficiently  rapid  on  the  time  scale  of  our  experiment, 

1 

then  intensity  contained  in  the  6sq  feature  will  lead  to  an  overestimated  lower  limit  for  binding 

energy,  consistent  with  the  experimental  results.  Such  dynamical  factors  have  been  observed  for 
the  benzene(Ar)j  cluster.45 

Figure  5  shows  the  benzyl  radical!  Ar)j  cluster  structure.  The  argon  atom  is  calculated  to  lie 

directly  above  the  aromatic  ring  system  and  slightly  toward  the  a-carbon  moiety.  This  structure  is 
similar  to  that  of  the  benzene! Ar)j  cluster.45 

A  further  test  of  our  empirical  model  is  to  account  for  the  van  der  Waals  vibrations  present 
in  the  benzyl  radical(Ar)j  cluster  spectrum.  Table  VI  tabulates  the  results  of  such  calculations. 

The  three  van  der  Waals  cluster  modes  are  determined  in  the  harmonic  limit  by  numerically 
mapping  out  the  one-dimensional  potential  energy  functions,  fitting  the  potentials  to  a  quadratic 
function,  and  calculating  the  reduced  mass  from  the  benzyl  radical  and  argon  entities.  In  the 
notation  of  Table  VI,  the  plane  of  the  benzyl  radical  is  the  xy-plane  with  the  y-axis  being  the 
principle  axis  of  the  benzyl  radical.  A  vibration  of  the  argon  atom  against  the  aromatic  ring  plane 
(z-axis  motion)  is  seen  to  yield  the  largest  vibrational  energy  mode  at  ca.  43  cm'1,  with  x-  and  y- 
translational  motion  parallel  to  the  aromatic  ring  plane  results  in  vibrations  at  11  and  10  cm'1, 
respectively.  The  bottom  portion  of  Table  VI  assigns  the  observed  vibrations.  Reasonable  overall 

agreement  exists  between  the  calculated  and  observed  van  der  Waals  modes. 

3 .  Benzyl  Radical fN2)  j 


The  MRES  of  the  benzyl  radical(N2)j  cluster  is  presented  in  Figure  6.  It  is  similar  to  that 
of  the  benzyl  radical(Ar)1  spectrum  with  a  few  important  and  interesting  differences:  the  relative 
intensities  for  the  A2  and  6aQ  transitions  are  smaller  than  for  the  bare  benzyl  radical  and  benzyl 

radical(Ar)j  cluster;  the  observed  spectroscopic  shifts  (see  Table  IV)  are  smaller  than  for  the 

benzyl  radicaKAr^  cluster,  intense  van  der  Waals  mode  structure  is  observed  at  the  1  A2  <—  1  B2 

origin  and  at  the  Aj  vibronic  feature;  and  this  van  der  Waals  mode  structure  is  different  for  the 

1  A2  and  Aj  transitions.  We  will  discuss  each  of  these  observations  in  turn  below. 
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The  decrease  in  intensity  of  the  A2  and  6aQ  features  relative  to  the  1  A2  Oq  and  A]  features 
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may  suggest  that  the  cluster  undergoes  dissociation  at  ca.  400  cm  above  the  1  A2  0  state.  One 
would  not  expect  that  such  an  intensity  decrease  would  result  from  changes  in  Franck-Condon 
factors  for  these  two  transitions,  considering  the  results  for  the  other  clusters  discussed  in  this 
report.  The  calculated  binding  energies  for  the  benzyl  radical(N2)j  cluster  are  given  in  Table  V. 

Since  the  calculated  binding  energy  for  the  cluster  is  ca.  530  cm'1,  the  potential  energy  function 

used  apparently  overestimates  the  binding  energy  by  roughly  100  cm'1. 

Table  V  lists  derived  spectroscopic  shifts  for  the  benzyl  radical(N2)j  cluster  system.  The 
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calculated  shift  for  the  1  A2  0q  origin  feature  of  -32.4  cm  is  a  factor  of  two  larger  than  the 

observed  shift  of  -16.6  cm'1  (Table  IV),  and  the  shift  for  A2  of  -5.78  cm'1  is  also  larger  than  the 

observed  shift  of  -0.7  cm'1  (Table  IV).  The  calculations  do  correctly  predict  the  trend,  however, 
yielding  a  smaller  shift  for  the  vibronically  mixed  states.  Table  V  shows  that  the  electrostatic  part 

of  the  intermolecular  potential  for  this  cluster  is  important,  contributing  a  total  of  20%  to  the 
binding  energy,  as  is  to  be  expected  for  N2  because  of  its  large  quadruple  moment.44 

The  calculated  1  B2  benzyl  radical(N2)j  cluster  structure  is  shown  in  Figure  7.  The 

nitrogen  molecule  is  calculated  to  lie  above  the  plane  of  the  benzyl  radical  with  its  bond  axis  parallel 
to  the  benzyl  radical  x-axis  (see  Table  III).  Calculation  of  the  structure  of  the  benzyl  radical(N2)j 
cluster  in  the  1  A2  electronic  state  shows  the  same  general  structure  as  the  ground  state  cluster 

except  that  the  nitrogen  molecule  is  displaced  by  0.07  A  along  the  positive  y-axis  (see  Table  III). 
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A  similar  calculation  for  the  cluster  structure  in  the  2  B2  electronic  states  yield  a  geometry  in  which 

the  nitrogen  molecule  bond  axis  is  parallel  to  the  benzyl  radical  y-axis  (see  Table  III  and  Figure  7). 

These  calculated  results  are  qualitatively  consistent  with  the  observed  van  der  Waals  vibrational 
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mode  structure  built  on  the  1  A2  Oq  and  Aj  features  of  Figure  6:  the  spectroscopic  structure 
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observed  for  the  1  A2  origin  corresponds  to  a  vibrational  van  der  Waals  mode  progression  and  the 
structure  at  A j  is  reminiscent  of  internal  rotation  such  as  that  observed  and  assigned  for  the 
benzene(N2)j  cluster.46  The  rotor  transitions  observed  here  may  also  be  associated  with  relatively 
free  rotation  of  the  N2  molecule  above  the  plane  of  the  benzyl  radical  as  we  now  discuss. 
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A  quantitative  interpretation  of  the  spectroscopic  structure  at  the  1  A2  electronic  origin  and 
Aj  features  of  Figure  6  can,  in  principle,  be  accomplished  through  the  use  of  empirical  potential 

calculations.  The  solid  lines  in  Figure  8  show  the  potential  energy  functions  generated  by 
translating  the  N2  molecule  along  the  ±  y-axis  relative  to  its  calculated  minimum  energy  location. 
The  small  displacement  between  the  minima  of  the  DQ(12B2)  and  Dj(12A2)  potentials  of  0.07  A  is 

barely  noticeable.  The  dashed  lines  correspond  to  the  two  quadratic  potentials  which  best  fit  the 
two  exponential-6-Coulomb  atom-atom  potentials  (Eq.  (1))  near  their  minima.  The  calculated 

potential  curves  are  generated  by  optimizing  all  but  the  one  degree  of  freedom  of  interest.  The 
fitted  functions  used  are  of  the  form  V=(l/2)k  Ay2,  in  which  k  for  both  the  12B2  and  12A2  states  is 

coincidentally  found  to  be  100  cm'VA2.  The  harmonic  vibrational  levels  using  this  value  for  k  and 
the  reduced  mass  for  this  translational  motion  [i.e.,  MBRMN2/(MBR  +  MN2)]  are  shown  in  Figure 

8.  These  derived  (harmonic)  potentials  are  employed  to  calculate  the  Franck-Condon  intensities 
displayed  in  Figure  6.  The  calculated  vibrational  spacing  of  12.5  cm*1  is  a  factor  of  two  lower 
than  the  observed  value  of  25.0  cm'1,  but  the  relative  transition  intensities  are  close  to  those 
observed. 

A  similar  calculation  can  be  performed  for  the  van  der  Waals  structure  built  on  the  Aj 

feature  (2  B2  mixed  vibronic  state)  shown  in  Figure  6.  Due  to  the  rotational  displacement  of  the 
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N2  molecule  in  the  2  B2  electronic  state  with  respect  the  ground  1  B2  state,  a  potential  energy 


function  corresponding  to  rotation  of  the  N2  molecule  is  generated.  The  solid  lines  in  Figure  9 
show  the  rotational  potential  energy  function  obtained  for  N2  rotation  in  the  1  B2  and  2  B2 

electronic  states.  A  two-fold  potential  is  seen  to  exist  with  barriers  of  13.0  and  17.2  cm'1, 
respectively.  The  molecular  symmetry  group  which  characterizes  this  two-fold  potential  is  the  G2 

group.  The  G2  group  is  isomorphic  to  the  C§  molecular  point  group:  details  of  the  generation  of 

this  molecular  symmetry  group  as  well  as  the  labeling  of  the  rotor  levels  are  given  elsewhere.47’48 
The  dashed  line  in  Figure  9  is  the  2  B2  potential  which  best  reproduces  the  observed  transitions. 

Table  VII  lists  the  observed  and  calculated  transitions  within  this  potential  model.  The  agreement 

is  remarkable,  if  not  fortuitous,  between  the  experimental  and  calculated  potential  functions. 

4.  Benzyl  Radical(CH4)j 

The  benzyl  radical(CH4)1  cluster  spectrum  is  presented  in  Figure  10.  The  intensity  of  the 
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1  A2  0q  and  6aQ  features  of  Figure  10  compared  to  those  of  the  bare  radical  (Figure  3),  as  well  as 

the  complete  lack  of  higher  energy  features  such  as  A3,  suggests  the  cluster  binding  energy  lies 
between  440  and  650  cm'1.  The  calculated  binding  energy  from  Table  V  of  ca.  540  cm"1  for  this 
cluster  is  consistent  with  the  experimental  results. 

The  observed  spectral  shifts  for  the  benzyl  radical(CH4)j  cluster  are  listed  in  Table  IV:  they 

are  roughly  twice  those  observed  for  the  argon  cluster,  and  three  times  those  of  the  nitrogen 
cluster.  The  calculated  shifts  given  in  Table  V  agree  well  with  the  experimentally  observed  values. 
The  vibronically  mixed  states  have  smaller  shifts  than  the  pure  1  A2  states  and  are  of  the  correct 

magnitude. 
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The  1  A2,  Aj,  A2,  and  6ag  features  of  Figure  10  all  contain  rich  spectroscopic  structure 

associated  with  low  energy  van  der  Waals  motion.  The  type  of  motion  giving  rise  to  such  structure 

is  likely  a  combination  of  methane  rotation  and  vibration,  such  as  that  observed  for  the 
benzene(CH4)j  cluster.  Empirical  potential  energy  calculations  for  benzene(CH4)j  show  that 

methane  is  situated  on  the  aromatic  ring  system  with  three  hydrogens  positioned  toward  carbon- 
carbon  bonds  giving  rise  to  a  cluster  of  effective  C3v  symmetry  49  Calculated  low  energy  van  der 


Waals  motion  for  this  cluster  consists  of  methane  rotation  described  by  a  three-fold  potential 
energy  function  plus  cluster  vibrational  motion.  Methane  contains  both  three-fold  and  two-fold 
rotational  axes.  Insight  into  the  preferred  type  of  rotor  motion  may  be  obtained  from  calculated 
cluster  equilibrium  structures  in  the  various  electronic  states.  Figure  1 1  shows  these  calculated 
structures.  Both  the  1  B2  and  l  Aj  electronic  states  have  similar  asymmetric  (double-well) 

structures  with  methane  located  near  the  a-carbon  moiety.  The  calculated  methane  configuration 
for  these  two  electronic  states  is  such  that  neither  a  pure  "two  hydrogen  down"  nor  "three 
hydrogen  down"  structure  exists.  Internal  rotational  motion  of  the  methane  molecule  will  most 

likely  occur  but  the  axis  of  its  rotation  is  uncertain.  The  double-well  potentials  for  methane  in  the 
benzyl  radical  1  B2  and  1  A2  states  may  greatly  effect  the  observed  spectra.  For  the  2  B2 

electronic  state  a  symmetric  (single-well)  "three  hydrogen  down"  methane  structure  is  calculated. 
The  expected  internal  rotational  motion  for  this  latter  state  would  be  one  arising  from  a  three-fold 
potential  energy  term. 
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Table  Villa  lists  the  observed  spectroscopic  features  for  the  1  Oq  and  6a0  regions  with 

assignments  based  on  either  a  three-fold  (G3  molecular  symmetry  group)  or  two-fold  (G2 

molecular  symmetry  group)  potential  energy  function.  The  top  portion  of  Table  Villa  presents 
results  for  the  benzyl  radical(CH4)j  cluster  spectrum  of  Figure  10  and  the  bottom  portion  of  the 
Table  Villa  lists  results  for  deuterated  methane.  Examining  first  the  G3  assignments,  four  out  of 

six  observed  features  can  be  assigned  as  methane  rotational  transitions;  the  remaining  two  features, 
by  default,  are  assigned  as  cluster  vibrational  modes.  Credibility  to  the  15.0  and  13.5  cm'1 

assignments  is  obtained  by  calculating  the  reduced  mass  isotope  correction  factor, 
[(Ir/M-dJ  ^  =  0-91,  which  correctly  predicts  the  isotopic  shift  for  a  cluster  vibrational  mode  treated 

as  an  uncoupled  harmonic  oscillator.  In  a  similar  fashion,  the  G2  assignment  of  the  same  features 
can  account  for  all  but  one  of  the  observed  features  as  arising  from  methane  rotation.  Both  the  G3 
and  G2  assignments  contain  only  the  barrier  height  as  an  adjustable  parameter.  Unfortunately,  the 

rotational  constant  for  methane  is  identical  for  both  three-fold  and  two-fold  rotor  motion. 


A  similar  ambiguity  in  feature  assignment  arises  for  motion  accompanying  the  Aj  and  A2 
features  of  Figure  10.  Table  VUIb  lists  the  observed  features  and  assignments.  Ironically,  the 
same  barrier  height  employed  in  the  calculations  presented  in  Table  Villa  fits  the  spectroscopic  data 

listed  in  Table  VMb.  Fewer  features  are  observed  for  these  vibronically  mixed  states,  nonetheless 
the  same  situation  arises:  both  G3  and  G2  assignments  can  be  made  of  the  data. 

The  empirical  potential  energy  calculations  can  be  used  in  a  straightforward  way  to  extract 
barrier  heights  for  three-fold  and  two-fold  rotor  motion.  The  calculations  should  be  able  to 
distinguish  which  type  of  motion  is  occurring.  The  lowest  barriers  found  for  three-fold  rotational 
motion  are  21.5,  32.6,  and  136.5  cm'1  for  the  12B2,  12A2,  and  22B2  electronic  states, 

respectively.  Two-fold  barriers  in  excess  of  200  cm'1  are  obtained  for  the  same  three  electronic 
states.  The  calculations  are,  therefore,  unable  to  predict  accurate  barrier  heights  for  the  rotor 
motion.  Two  important  results  can  be  extracted  from  the  above  discussion.  First,  the  calculations 
succeed  in  predicting  the  benzyl  radical(N2)1  rotational  barrier  is  not  repeated  for  the  methane 

cluster  system.  Second,  and  equally  important,  care  must  be  exercised  in  extracting  a  rotational 
barrier  from  spectroscopic  data  which  include  low  energy  van  der  Waals  vibrational  modes. 

5 .  Benzyl  Radical  (C2H6)  j 

The  benzyl  radical  (C2H6)j  cluster  spectrum  is  presented  in  Figure  12.  The  A}  feature  is 

seen  to  be  much  more  intense  than  the  remaining  features.  In  addition,  the  presence  of  the  A3 

feature  suggests  a  binding  energy  greater  than  650  cm'1  exists  for  this  cluster.  This  is  consistent 

with  the  calculated  binding  energies  listed  in  Table  V.  Two  stable  conformations  are  calculated  and 

are  pictorially  represented  in  Figure  13.  The  more  stable  of  the  two  conformers  of  this  cluster, 

labeled  1,  has  ethane's  carbon-carbon  bond  axis  nearly  parallel  to  the  benzyl  radical  principle  axis. 

Structure  2  has  the  ethane  carbon-carbon  bond  nearly  parallel  to  the  benzyl  radical  x-axis  (see 

Table  III).  In  Figure  12  we  have  assigned  the  well-resolved  spectroscopic  structure  to  the  lowest 
energy  conformer  1.  The  broad  feature  in  the  vicinity  of  the  Aj,  A2,  and  6aQ  transitions  be  may 


due  to  van  der  Waals  vibrational  structure  or  the  presence  of  transitions  associated  with  the  cluster 
structures  labeled  2. 

The  spectroscopic  shifts  corresponding  to  conformation  1  are  listed  in  Table  IV.  The 
corresponding  calculated  shifts  in  Table  V  overestimate  the  observed  shifts  by  roughly  10  cm'1  for 

o  1 

the  1  A2  states,  and  20  cm  for  the  vibronically  mixed  states. 

Low  energy  van  der  Waals  modes  are  observed  for  this  cluster  species  as  in  the  cluster 
systems  already  discussed.  Due  to  the  increased  complexity  of  the  spectrum  and  possible  existence 
of  multiple  conformations,  an  analysis  of  the  vibrational  van  der  Waals  modes  is  not  attempted. 

6.  Benzyl  radical(C3Hg)j 

The  benzyl  radical(C3Hg)j  cluster  spectrum  is  presented  in  Figure  14.  It  is  similar  to  that 
found  for  the  ethane  cluster  spectrum  of  (Figure  12)  in  that  the  feature  has  the  largest  intensity 
and  the  A3  feature  is  observed.  Additional  spectroscopic  structure  to  the  blue  of  the  A3  feature  up 

to  roughly  1000  cm"1  is  consistent  with  the  calculated  binding  energy  of  ca.  970  cm'1  listed  in 
Table  V.  No  features  are  observed  to  the  blue  of  the  spectrum  shown.  Three  stable  structural 
conformations  are  calculated  for  this  cluster,  with  ground  state  binding  energies  given  in  Table  V 
and  pictorial  representations  shown  in  Figure  15.  The  lowest  energy  structure,  1,  has  propane's 
three  carbon  atoms  lying  in  the  yz  plane  perpendicular  to  the  xy  plane  of  the  benzyl  radical  (see 
Table  ELI);  structure  2  has  propane's  three  carbon  atoms  lying  in  a  plane  parallel  to  the  plane  of  the 
benzyl  radical  (xy)  above  carbons  2,  3, 4,  and  5  of  the  benzyl  radical  (see  Table  III);  and  structure 
3  has  propane's  three  carbon  atoms  also  in  a  plane  parallel  to  that  of  the  benzyl  radical  xy  plane 
above  carbons  2  and  3  (see  Table  III).  Similar  results  have  been  found  for  the  toluene/propane 

cluster  system.50  A  structural  conformation  assignment  of  the  features  in  Figure  14  is  not  unique. 
A  plausible  scheme  can  label  the  most  intense  features  (1  A2,  Aj,  A2, 6a,  and  A3)  as  belonging  to 

structure  1,  and  the  weaker,  red-shifted,  features  to  a  second  conformation  such  as  2  (or  3). 
Such  difficulties  in  conformation  assignment  also  make  van  der  Waals  mode  assignment 
prohibitively  complex. 


V. 


CONCLUSIONS 


In  this  report  we  have  presented  spectra  and  calculations  for  the  benzyl  radical  solvated 
with  the  non-polar  solvents  Ar,  N2,  CH4,  C2H6,  and  C3Hg.  The  experimentally  observed 

solvation  effects  of  the  benzyl  radical  are  not  dissimilar  from  those  found  for  the  benzene 
system :45.46,49  smajj  SpeCQ-ai  re(j  shifts,  the  presence  of  van  der  Waals  modes,  and  binding 
energies  in  the  range  of  400  to  950  cm'1. 

The  varied  solvents  studied  here  allow  a  test  of  empirical  atom-atom  based  potential  energy 
calculations  to  be  performed.  Binding  energies  for  the  argon  and  nitrogen  clusters  are 
overestimated  by  roughly  50-100  cm'1,  whereas  for  the  hydrocarbon  solvents  the  calculated 
binding  energies  are  consistent  with  the  experimental  results.  The  predicted  trends  of  the 
spectroscopic  shifts  for  all  the  solvents  examined  are  accurately  modeled,  but  the  size  of  the  shifts 
are  incorrect  by  as  much  as  a  factor  of  three.  The  calculated  van  der  Waals  low  energy 
vibrational/torsional  modes  yield  mixed  results:  for  argon  and  nitrogen  the  calculations  are  of 
predictive  value  while  for  methane  the  derived  internal  rotational  barriers  are  too  large.  Certainly 
there  is  room  for  improvement  as  far  as  atom  based  van  der  Waals  parameters  are  concerned.  An 
interesting  area  to  explore  theoretically  is  the  possibility  of  ab  initio  derived  atom-based  van  der 
Waals  parameters.  If  such  a  theoretical  advance  is  made,  better  modeling  of  excited  states  (without 
assumptions  about  atomic  van  der  Waals  parameters)  would  become  tractable. 
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Table  I.  Non-bonded  Interaction  Parameters  for  the  DQ  (12B2),  Dj  (12A2),  and  D2 
(22B2)  Electronic  States. 


Species 

A  (kcal/mol) 

B  (A'1) 

C  (A6  kcal/mol) 

Benzyl  Radical  DQ 
/  Solvent  Parmeters: 

Carbon3 

30026. 

3.3845 

472.2 

Hydrogen3 

1687.6 

3.537 

43.42 

Nitrogen*5 

60862. 

3.7801 

329.5 

Argonc 

40103. 

3.1258 

1577. 

Benzyl  Radical  Dj 
Parameters:*1 

C 

34620. 

3.3845 

544.5 

H 

1687.6 

3.537 

43.42 

Benzyl  Radical  D2 

Parameters. 

C 

32578. 

3.3845 

512.3 

H 

1687.6 

3.537 

43.42 

a  Ref.  40. 
b  Ref.  43. 
c  Ref.  51. 

d  These  parameters  are  obtained  by  adjusting  the  carbon  parameters  so  as  to  reproduce  the 
experimentally  measured  benzyl  radical(Ar)j  spectroscopic  shifts.  See  text  and  Tables  IV 

and  V. 


Table  II.  Ab  Initio  Calculations  of  the  Excited  States  of  the  Benzyl  Radical. 


Transition  Energy  Estimation: 

State  Symmetry 

Root3 

Energy  (a.u.) 

RHF  optimized  geometry13  CASSCF  optimized 

geometry13 

B2 

first 

-269.128883 

(0.0) 

-269.128961 

(0.0) 

B2 

second 

-269.017838 
(24374  cm'1) 

-269.024851 
(22852  cm'1) 

*2 

first 

-269.012503 
(25545  cm'1) 

-269.021319 
(23627  cm'1) 

3  Active  space  consists  of  7  tc  orbitals  (Dunning-Huzinaga  [9s,5p]  double  zeta  valence  basis  set). 
RHF  result:  (b2)2  (b2)2  (a2)2  (b2)1  (a^  (b2)  (b2),  B2  state  symmetries  resulted  in  404 
configurations,  and  A2  state  symmetries  resulted  in  380  configurations. 

b  Energy  gradients  reduced  to  less  than  10"4  hartree/bohr. 


Ionization  Energy  Estimation: 

UHF  B2  ground  state  energy  =  -269.071 160 
RHF  Aj  +  1  cation  energy  =  -268.801410 

Calculated  I.P.  =  59,210  cm'1,  observed  I.P.  =  58,358  cm'1  (See  Ref.  9). 


Table  III.  Calculated  Partial  Charges  for  the  Benzyl  Radical  in  the  12B2>  12A2,  and  22B2  Electronic 
States.3  The  number  of  grid  points  used  for  each  electronic  state  is  5 19. 


Benzyl  Radical  Axis  System13 
and  Numbering  Scheme: 


Electronic  State 

Atom 

Partial  Charge  (e') 

Calc.  Residual 

Charge  (e‘) 

1 

-0.553 

2 

0.243 

12B2  (D0) 

3 

-0.149 

0.000410 

4 

-0.222 

5 

-0.104 

6 

0.203 

7 

0.139 

8 

0.163 

9 

0.144 

1 

-0.547 

2 

0.063 

12Aj  (D,) 

3 

-0.023 

0.000917 

4 

-0.184 

5 

-0.167 

6 

0.207 

7 

0.107 

8 

0.144 

9 

0.150 

1 

"  -0.654 

2 

0.540 

22B2  (D2) 

3 

-C.350 

0.000744 

4 

-0.164 

5 

-0.041 

6 

0.214 

7 

0.175 

8 

0.150 

9 

0.105 

3  CASSCF  12B2 

optimized  geometry  used,  grid  spacing  =  2.0  bohr,  and  shell  thickness  = 

b  The  origin  of  the  axis  system  corresponds  to  the  center-of-mass  of  the  benzyl  radical. 

0 

Table  V.  Calculated  Binding  Energies  (Vgg  cm’*)  for  Benzyl  Radical  Clusters  in  Different 

Electronic  States  and  Derived  Spectroscopic  Shifts  for  the  Two  Observed  Transitions. 

v 

VBE 

—  V  + 

v  VDW 

Vc 

Cluster  Shift, 

AE 

12B2: 

Ar 

-399.5 

-399.5 

0 

N2 

-519.5 

-422.7 

-96.8 

ch4 

-509.2 

-493.3 

-15.9 

c2h6  (1) 

-708.7 

-707.9 

-0.8 

C2H6  <2) 

-666.6 

-685.6 

19.0 

C3H8  (1) 

-930.0 

-898.8 

-31.2 

C3H8  (2) 

-865.1 

-852.1 

-13.0 

C3H8  (3) 

-855.2 

-879.7 

24.5 

l2^:3 

Ar 

-423.3 

-423.3 

0 

-23.8 

N2 

-551.9 

-447.9 

-103.9 

-32.4 

ch4 

-544.3 

-519.5 

-24.8 

-35.1 

C2H6b 

-748.2 

-753.3 

5.1 

-39.5 

C3H8b 

-976.3 

-957.1 

-19.2 

-46.3 

22B  2:a 

Ar 

-412.9 

-412.9 

0 

-13.4 

N2 

-525.2 

-436.0 

-89.2 

-5.7 

ch4 

-534.3 

-540.7 

6.4 

-25.1 

C2H6b 

-751.1 

-723.4 

-27.7 

-42.4 

C3H8b 

-971.2 

-920.4 

-50.8 

-41.2 

0  o 

aThe  binding  energies  for  the  excited  states  are  calculated  by  adjusting  the  1  A2  and  2  B2 
carbon  parameters  to  reproduce  spectroscopic  shifts  for  benzyl  radical(Ar)j.  The  calculated 

shifts  for  the  remaining  clusters  are  obtained  using  the  modified  excited  state  van  der  Waals 
parameters  and  appropriate  partial  charges  for  each  electronic  state  (see  Table  I). 

b  Only  lowest  energy  conformer  result  is  presented  here. 

Table  VI.  Benzyl  Radical(Ar)1  Calculated  van  der  Waals  Vibrations  and  Assignments. 


Mode  Calculation: 

mode  (cm'1) 

Electronic  State 

rx 

ly 

<z 

1^2 

11.0 

10.1 

43.4 

i2a2 

11.4 

10.4 

43.9 

2\ 

11.2 

10.4 

43.4 

vdW  Mode  Assignment: 

Observed 

Calculated5 

i2a2 

41.2 

lz 

43.9 

A, 

6.9 

ty 

10.4 

22.4 

lx  ly 

V 

21.8 

38.9 

41.6 

43.4 

‘I 

43.4 

^2 

6.0 

ly 

10.4 

14.5 

11.2 

6a1 

13.0 

tx 

11.4 

40.2 

t2 

43.9 

a  Plane  of  aromatic  ring  is  the  xy-plane,  the  y-axis  is  the  principle  axis  (see  Table  III). 

5  Calculated  modes  for  the  12A^  and  6a1  states  are  obtained  using  ^A^  state  van  der  Waals 

parameters,  whereas  Aj  and  A^  state  modes  are  obtained  using  22B2  state  van  der  Waals 
parameters. 


Table  VII.  Benzyl  Radical(N2)j  Internal  Rotational  Transitions. 


Observed  Transitions3 

Assignment 

Calculated  Transitions 

0.0 

*C$ 

O 

i 

O 

0.0 

1.0 

T 

o 

0.2 

6.0 

2a’  -e-  0a' 

8.6 

11.0 

3a"  0a’ 

10.8 

15.0 

4a'  <-  0a' 

14.4 

17.0 

vdw  vib. 

22.0 

5a"  <-  0a‘ 

22.0 

24.0 

6a'  <-  0a’ 

22.6 

3  Transitions  are  relative  to  the  Aj  feature  located  at  22,322.5  cm'1. 
b  The  potential  function  is  V=V2/2(l-cos2<t>)  with  V2=18.0  cm"1,  B=2.0  cm 


o  0  1 

Table  Villa.  Benzyl  Radical  (methane^  Low  Energy  Modes  for  the  1  A2  Oq  and  6aQ  Transitions. 


Benzyl  Radical(CH4)j 


Observed  (cm'1)2 

G3  assignment15 

Calc. 

G2  assignment0 

Calc. 

0.0 

0a  4-  0a 

0.0 

0a*  4-  0a' 

0.0 

4.8 

le  4-  le 

5.2 

la”  4-  0a' 

4.7 

15.0 

vdw.  vib. 

- 

2a'  4-  la" 

15.1 

19.8 

2e  4—  le 

20.9 

2a’  4-  0a' 

20.8 

38.7 

vdw.  vib. 

- 

3a"  4-  la" 

41.8 

45.9 

3a  4-  0a 

46.8 

3a"  4-  0a’ 

46.8 

Benzyl  Radical(CD4)1 


Observed  (cm'1)d 

G3  assignment15 

Calc. 

G2  assignment0 

Calc. 

0.0 

0a  4-  0a 

0.0 

0a'  4-  0a' 

0.0 

5.0 

le  4-  le 

2.6 

2a'  4-  la" 

7.3 

9.9 

2e  4-  le 

10.5 

2a'  4-  0a' 

10.5 

13.5 

vdw.  vib . 

- 

vdw.  vib. 

- 

20.1 

3a  4—  0a 

23.4 

3a"  4-  la" 

20.4 

24.3 

vdw.  vib. 

- 

3a"  4-  0a' 

23.5 

a  Relative  to  12A2  electronic  origin  located  at  21953.7  cm'1. 
b  Calculated  using  V3/2(l-cos3<]))  potential  with  V3=3.0  cm'1  and 
B=5.2  cm'1  for  CH4  or  B=2.6  cm'1  for  CD4. 
c  Calculated  using  V2/2(l-cos2<}>)  potential  with  V2=2.0  cm'1  and 
B=5.2  cm'1  for  CH4  or  B=2.6  cm'1  for  CD4. 
d  Relative  to  l2^  electronic  origin  located  at  21954.5  cm'1. 


Table  Vlllb.  Benzyl  Radical  (methane)  j  Low  Energy  Modes  for  the  Aj  and  A2  Electronic 
Transitions. 


Benzyl  Radical(CH4)j 


Observed  (cm'1)3 

G3  assignment15 

Calc. 

G2  assignment0 

Calc. 

0.0 

0a  <—  0a 

0.0 

0a'  0a’ 

0.0 

4.0 

le  <—  le 

5.2 

la"  <-  0a‘ 

4.7 

13.9 

vdw.  vib. 

- 

2a’  la" 

15.1 

18.9 

2e  «-  le 

20.9 

2a’  0a' 

20.8 

Benzyl  Radical(CD4)1 


Observed  (cm'1)11 

G3  assignment15 

Calc. 

G2  assignment0 

Calc. 

0.0 

0a  <—  0a 

0.0 

0a'  <-  0a' 

0.0 

4.0 

le  <—  le 

2.6 

2a'  <-  la" 

7.3 

9.1 

2e  <—  le 

10.5 

2a' <-  0a' 

10.5 

20.1 

3a  <—  0a 

23.4 

3a"  <-  la" 

20.4 

3  Relative  to  Aj  feature  located  at  22294.0  cm’1. 
b  Calculated  using  V3/2(l-cos3<j))  potential  with  V3=3.0  cm'1  and 
B=5.2  cm’1  for  CH4  or  B=2.6  cm’1  for  CD4. 
c  Calculated  using  V2/2(l-cos2(())  potential  with  V2=2.0  cm'1  and 
B=5.2  cm*1  for  CH4  or  B=2.6  cm'1  for  CD4. 
d  Relative  to  Aj  feature  located  at  22293.9  cm’1. 


FIGURE  CAPTIONS 


Figure  1.  Optical  arrangement  for  MRES  study  of  radicals. 

Figure  2.  Electronics  arrangement  for  MRES  study  of  radicals. 

Figure  3.  The  benzyl  radical  MRES.  The  features  labeled  Aj  and  A3 

arise  from  vibronic  coupling  between  the  non-totally  symmetric  6b  vibration  of  the  1 
electronic  state  and  the  2  B2  electronic  origin.  The  feature  labeled  is  assigned  as  the 

1  9 

18b  vibration  of  the  1  A2  electronic  state  which  vibronically  couples  and  borrows 
intensity  from  the  2  B2  electronic  origin.  (See  reference  19) 

Figure  4.  The  benzyl  radical(Ar)^  spectrum.  Note  that  the  spectrum  ends  at  roughly  450 

1  9  0 
cm  above  the  1  A^  0q  transition. 

Figure  5.  Calculated  structure  of  the  benzyl  radical(Ar)j  cluster  in  the  DQ  1  B2 
electronic  state.  The  argon  atom  lies  above  the  plane  at  z  =3.50  A  and  y  =  -.19  A.  The 
circle  representing  the  argon  atom  corresponds  to  the  argon  van  der  Waals  radius  (~1 .8 

A). 

Figure  6.  The  benzyl  radical(N2)1  spectrum.  The  calculated  spectrum  corresponds  to  a 
vibration  of  the  N2  molecule  along  the  plane  of  the  benzyl  radical.  See  text  for  details  of 
the  calculation. 

Figure  7.  Calculated  structure  of  the  benzyl  radical (N2)j  cluster  in  the  DQ  12B2  and  D2 
2 

2  B2  electronic  states.  The  bond  axis  of  the  nitrogen  molecule  is  parallel  to  the  benzyl 
radical  x-axis  in  DQ  and  perpendicular  to  it  in  D2  (see  Table  ID). 

Figure  8.  The  calculated  potential  energy  surfaces  corresponding  to  a  van  der  Waals 
vibration  of  the  N2  molecule  along  the  plane  of  the  benzyl  radical  are  shown  for  the  DQ  and 
Dj  electronic  states.  Specifically,  the  one-dimensional  functions  represent  the 
displacement,  from  its  minimum  energy  location,  of  the  N2  entity  along  the  y-axis  shown 


in  Table  HI. 


Figure  9.  The  calculated  potential  energy  surfaces  corresponding  to  the  N2  rotation 
above  the  plane  of  the  benzyl  radical  in  the  DQ  and  D2  electronic  states  (solid  lines). 

The  dashed  line  represents  the  fitted  potential  which  best  reproduces  the  spectroscopic 
data.  See  Table  VII  for  the  explicit  form  of  the  potential  used. 

Figure  10.  The  benzyl  radical(CH4)1  MRES.  Note  that  the  spectrum  ends  at  ca.  450 

1  9  0 

cm  above  the  1  Aj  0q  transition. 

Figure  11.  Pictorial  representation  of  the  benzyl  radical(CH4)j  cluster  structure  in  the 
Dq,  Dj,  D2  electronic  states  of  the  radical.  Note  the  change  in  methane  orientation  between 
the  1  2B2  and  1  2A2  states,  and  the  2  2B2  state. 

Figure  12.  The  benzyl  radical(C2H6)j  MRES.  Note  that  the  spectrum  ends  at  ca.  650 

-i  i  o 

cm'  above  the  1  Oq  transition. 

Figure  13.  Two  structures  for  the  benzyl  radicaKCjH^j  in  the  DQ  electronic  state.  The 

C-C  bond  axis  of  the  ethane  is  nearly  parallel  to  the  y-  or  x-axis  of  the  benzyl  radical  for 
Structures  1  and  2,  respectively. 

Figure  14.  The  benzyl  radical(C3Hg)j  MRES  .  Note  that  the  spectrum  extends  beyond 

1  9  0 

900  cm'1  above  the  1  A ^  0q  transition. 

Figure  15.  Three  structures  for  the  benzyl(C3Hg)j  cluster  in  the  DQ  electronic  state. 
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